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Abstract 

We finely describe the speed of "coming down from infinity" for birth and death 
processes which eventually become extinct. Under general assumptions on the birth and 
death rates, we firstly determine the behavior of the successive hitting times of large 
integers. We put in light two different regimes depending on whether the mean time for 
the process to go from n + 1 to n is negligible or not compared to the mean time to reach 
n from infinity. In the first regime, the coming down from infinity is very fast and the 
convergence is weak. In the second regime, the coming down from infinity is gradual and 
a law of large numbers and a central limit theorem for the hitting times sequence hold. By 
an inversion procedure, we deduce that the process is a.s. equivalent to a non-increasing 
function when the time goes to zero. Our results are illustrated by several examples 
including applications to population dynamics and population genetics. The particular 
case where the death rate varies regularly is studied in details. 

Key words: Birth and death processes. Coming down from infinity. Hitting times. Central 
limit theorem. 

MSC 2010: 60J27, 60J75, 60F15, 60F05, 60F10, 92D25. 

1 Introduction and main results 

Our goal in this paper is to finely describe the "coming down from infinity" for a birth 
and death process. We are motivated by the study of population dynamics and population 
genetics models with initially large populations. For this purpose, we first decompose the 
trajectory of the process with respect to the hitting times of large integers. We then study 
the small time behavior of the continuous time process when it comes down from infinity. 

The population size is modeled by a birth and death process > 0) whose birth 

rate (resp. death rate) at state n G N is (resp. /U„). In the whole paper, the rates are 
nonnegative and the rates Hn are positive for n > 1. Moreover, we assume that /tq = Aq = 0 
for practical purpose. The latter implies that 0 is an absorbing state. Such processes have 
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been extensively stndied from the pioneering works on extinction [9] and qnasi-stationary 
distribntion [19]. 


It is well known [9, 10] that 


V — 

* ^ \ in 


i>\ 




= oo 


( 1 ) 


is a necessary and snfhcient condition for almost snre absorption of the process at 0, where 


TTi = — and for n >2, 
fJ-i 


1 


Under Condition (1), we hrst dehne the law Pqo of the process starting from inhnity with 
valnes in N U {oo} (see Lemma 2.1) as the limit of the laws of the process issned from n. 

When the limiting process is non-degenerate, it hits hnite valnes in hnite time with positive 
probability. This behavior is captnred by the notion of "coming down from inhnity". A key 
role is played by the decreasing seqnence (T„)n>o of hitting times dehned as 


Tn := infjt > 0, A(t) = nj. 


As proved in [17, p.384] and in [2, Chap.3], 

Ei(ro) = ^7ri and E„+i(r„) =-^ X] for n > 1. (2) 

Remark that in case of pnre-death process, the law of Tn nnder En+i is exponential with 
parameter //n+i and for n G N, E„_|_i(T„) = 

Characterizations of the coming down from inhnity have been given in [3, 7]. They rely on 
the convergence of the mean time of absorption when the initial condition goes to inhnity or 
eqnivalently to the convergence of the series 


S= lim E„(T„) = ^,r. + j;— ^ = E E 


A 


n+1 


• • Aj i 


i>l 


Atj T! 
n>l i>n+l 


n>0 \i>n+l 


/^n+1 ' ' ' 


< -1-00. 


(3) 


This is eqnivalent to the existence and nniqneness of the qnasi-stationary distribntion at 0 (see 
[19], [7]) and to the hniteness of some exponential moments of Tq. Fnrthermore, monotonicity 
properties allow ns to show that this is also eqnivalent to instantaneons almost-snre coming 
down from inhnity (Proposition 2.5). 

In the whole paper, we snppose that Assnmption (1) holds and from Section 3 onward, we 
assnme that (3) is satished, that is, the process instantaneonsly comes down from inhnity. It 
gnarantees the hniteness of all moments of nnder En-i-i and nnder Eqo) for which we have 
an explicit expression (Proposition 2.2). In Section 3, we pnt in light two different regimes 
for the asymptotic behavior of T„/Eoo(7n), depending on whether the mean time to go from 
n -|- 1 to n is negligible or not compared to the mean time to reach n from oo. In the hrst 
regime, the coming down from inhnity is very fast and the limit is random. In the second one, 
the coming down is gradnal and dne to the accnmnlation of small independent contribntions, 
which leads to a law of large nnmbers. More precisely, we assnme that 

®'n-|-l(^n) _ XnTTn Si>n+1 ^ ^ 

lEoo(En) ^*>1+1 
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In the first (fast) regime (a > 0) and under the additional assumption, which is stronger than 

( 1 ): 

I := lim — < 1, (4) 

n^+oo 

we prove that T„/Eoo(7n) converges in law to a non-degenerate random variable whose dis¬ 
tribution is characterized by I and a. 

In the second (gradual) regime (a = 0), we prove a weak law of large numbers under the 
following second moment assumption 


sup 

n>0 


^n+ljTn) 

{En+l{TnW 


< -1-00. 


(5) 


More precisely, we prove that the sequence {Tn/Eoo{Tn)) converges in probability to 1. Under 
some additional variance assumptions, we also obtain a central limit theorem. 

Thanks to (2) and to forthcoming (8), we note that both expectations in (5) can be written in 
terms of the birth and death rates. Condition (5) is fulhlled in many cases we have in mind. 
For instance, it holds for pure death processes. 

We will see in the next section some more tractable conditions ensuring (1), (3), (4) and (5). 
In the second regime, under (5) and the following additional condition 



En+l{Tn) 

Eoo(Tn) 


2 


< + 00 , 


( 6 ) 


which means that the convergence of E„_|_i(T„)/Eoo(Tn) to 0 is fast enough, one also get a 
strong law of large numbers for Tn/Eoo (Tn). 


We then derive in Section 4 the small time behavior of the process X. We prove that 

t^o v{t) 

where v is the generalized inverse function of n i—>• Eoo(T'„) = Ylj>n ■ 

v{t) = inf{n > 0; Eoo(T„) < t}. 


The limit holds in probability in the hrst regime. Remark that in this fast case, Tn/Eoo{Tn) 
converges in law (an not in probability) to a random variable but nevertheless, X{t) behaves as 
v(t) for t small. That is due to the fact that Eoo(T'[„a;]) is negligible with respect to Eoo(Tn) for 
any x > 1 and for large n. In the second regime, one needs some additional assumptions and 
almost sure convergence can be obtained. The proof relies on two ingredients: the short time 
behavior of the non-increasing process equal to n on [r„,T„_i[ and the control of the height 
of the excursion of the process X during the time interval [r„,T„_i). Technical assumptions 
are required in the second regime to estimate the variations of Eoo(7n) and to deduce the 
behavior of X from that of {Tn)n, by a non trivial inversion procedure. Our motivations and 
applications from population dynamics and population genetics meet these assumptions. Thus, 
our results cover general birth and death models including many different ecological scenarios. 
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as competition models with polynomial death rates [16] or Allee effect [12], Lambert [13] 
characterizes the distribution of the absorption time for the logistic branching process starting 
from infinity. Our work extends in different way the case of Kingman coalescent, for which 
speed of coming down from infinity has already been obtained by Aldous [1]. More generally, 
in the gradual regime, the behavior of the process coming down from in finity is similar to that 
of A coalescent obtained in [4, 14], We also note that our approach relies on the decomposition 
of the trajectory of X with respect to the reaching times of the successive integers and our 
results could be extended to processes with several births. Another motivation for the results 
below is the study of birth and death process in random environment and in particular the 
study of the regulation of the population during unfavorable periods. This latter is a work in 
progress. 

The paper is organized as follows. In the next section, we work under the absorption assump¬ 
tion (1) and prove the existence of the law of the process starting from infinity. Thus we gather 
general characterizations of the coming down from infinity and we show the equivalence to the 
a.s. instantaneous coming down. Focusing in Section 3 on birth and death processes satisfying 
(4) or (5), we describe the hitting times of large integers. In Section 4, we obtain a law of large 
numbers describing the small time behavior of the process X. Examples and applications are 
provided all along the paper and illustrate the different regimes. The last Section 5 focuses 
on regularly varing death rate and provides our main application to small time behavior for 
population dynamics and population genetics one-dimensional processes coming down from 
infinity. 

2 Preliminaries and coming down from infinity 

2.1 Preliminaries 

The first lemma allows us to define the law of the process starting from infinity. It is based 
on monotonicity arguments following Donnelly [8]. We set N := {0,1,. ..} U {oo} and for any 
T > 0, we denote by D^([0, T]) the Skorohod space of cadlag functions on [0,T] with values 
in N. 

Lemma 2.1. Under (1), the sequence (Pn)n converges weakly in the space of probability mea¬ 
sures on D^([0, T]) to a probability measure Poo- 

At this point, the limiting process is not assumed to be finite for positive times. 

Proof. We follow the tightness argument given in the first part of the proof of Theorem 1 by 
Donnelly in [8]. Indeed, no integer is an instantaneous state for the process {Xn,Hn < oo for 
each n > 0) and the process is stochastically monotone with respect to the initial condition. 
It ensures that Assumption (Al) of [8] holds. In addition. Assumption (1) ensures that the 
process almost surely does not explode and (A2) of [8, Thm. 1] is also satisfied by denoting 
the birth and death process X issued from n and stopped in N. 

Then the tightness holds and we identify the finite marginal distributions by noticing that for 
A: > 1, for ti,..., > 0 and for oi,..., G N, the quantities P„(Ai(ti) < ai, • • • ,X{tk) < a^) 

are non-increasing with respect to n G N (and thus converge). □ 
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We now focus on the time spent by the process > 0) to go from level n + 1 to level n. 

For n > 0, we introduce the function 

Gnia) := E„+i(exp(-aT„)), a > 0. 


Proposition 2.2. Suppose that (1) holds. For any a > 0 and n > 1, we have 

fin + a fin 1 


Gn{o) — 1 + 


^nGn—l{(l) 


Moreover, for every n > 0, 
2 


E„+i(T2) = , E„+i(r^) = —^ ^ Ai7rjEj+i(ri)Vari+i(ri). 


^n'^n 

i>n 


(7) 


( 8 ) 


Proof. We denote by Tn a random variable distributed as Tn under E^+i and consider the 
Laplace transform of Tn- Following [3, p. 264] and by the Markov property, we have 


Tn—l — l{y„=—IjE'n + \y„=1} {E/u + 'Ll + Tn—l) 


where Yn, En, t^-i and are independent random variables, En is an exponential random 
variable with parameter + fin and is distributed as t„_i and E(l^ = 1) = 1 — E(l^ = 
-1) = A„/(An + fin)- Hence, we get 


i(n) — 


Y fin 

n 4“ A 72 Y fin 


Gn{a)Gn-i{a) 


A. 


+ 


fin 


An Y fin An + fin 


and (7) follows. 

Differentiating (7) twice at a = 0, we get 


En(r2_l) = —En+l(r2) + 2En(rn_l)2, n > 1. 
fin 

In the particular case when Aat = 0 for some N > n, a simple induction gives 

En+i(r 2 ) = Y^ ^ XmEi+i{Tif 


(9) 


and (8) is proved. In the general case, let N > n. Thanks to Assumption (1), Tq is finite 
and the process a.s. does not explode in finite time for any initial condition. Then Tn 
is finite and T^r —)> +oo En+i-a.s., where we use the convention = +oo on the event 
{Vt >0:X{t)^ A^}. The monotone convergence theorem yields 


En+l (T„;Tn < Ttv) > En+l(T^). 

N^+oo 

Let us consider a birth and death process X^ with birth and death rates (A^,p^ : k > 0) 
such that {X^,fi^) = {\k, hk) for A: / iV and = 0, fi^ = fiN- 
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Since {Xt : t < T]\f) and : t < T^) have the same distribution under P^+i) we get 
E„+i (r2;T„ < Tn) = En+i < T^) , 

which yields 

E„+i(t 2) = hm E „+1 ((r„^)2;r„^ < T^) < lim E„+i {{T^f) , 

N^oo N^oo 

where the convergence of the last term is due to the stochastic monotonicity of with 
respect to N under En+i- Since is stochastically smaller than Tn under En+i, we have 
also 

En+l{iTnf)>En+l{{T^f). 

We deduce that 

En+l{{Tnf) = lim En+l{{T^f) = hm ^ ^ 

N^oo N^oo A„7r„ 

where the last identity comes from (9). Adding that Ej+i(T^^) is non-decreasing with respect 
to N yields the expected expression for Eri+i((rn)^) by monotone convergence. 

The third moment is obtained similarly by differentiating (7) three times, which gives the 
recurrence equation 

EniT^_i) = —En+i{T^) + 6E„(T„_i)Var„(r„_i), n > 1. 

fJ-n 

The coupling argument we have used above allows us to conclude. □ 


Remark 2.3. Using Proposition 2.2, (5) writes 


AjVTj / Ej_|_i(Tj) 
n>0 ^n'^n \E)i-|-l(T^) 


2 

< - 1 - 00 . 


( 10 ) 


2.2 Instantaneous coming down from infinity 

We now dehne a strong notion of coming down from inhnity corresponding to the behavior of 
birth and death processes under (1) and (3): the process comes down instantaneously almost 
surely. 

Definition 2.4. The process {X(t),t > 0) instantaneously comes down from infinity if for 
any t > 0, 

lim lim < t) = 1. (11) 

m—^oo k^+oo 

Note that (11) is equivalent to 

Eoo(Vt > 0,X{t) < -t-oo) = 1. 

Let us now show that (11) is satished under (1) and (3). In fact we give several necessary 
and sufficient conditions for {X{t),t > 0) to come down from inhnity. The hrst two ones 
are directly taken from [7]. We add here an exponential moment criterion. We also mention 
that it is equivalent to the existence (cf. [19]) and uniqueness (cf. [7]) of a quasi-stationary 
distribution for the process X. 
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Proposition 2.5. Under condition (1), the following assertions are equivalent: 

(i) The process {X(t),t > 0) instantaneously comes down from infinity. 

(ii) Assumption (3) is satisfied: S < +oo. 

(iii) supfc>oEfc[ro] < +CX). 

(iv) For all a > 0, there exists fea G N such that (exp(aTfc^)) < +cx). 

Proposition 2.5 implies in particular that under (1) and (3), the moments of Tn under and 
Poo are finite. Moreover their explicit expression can be derived from Proposition 2.2 and will 
be useful in the rest of the paper. 

Remark 2.6. It is proved in the next Lemma 2.1 (i), that the equivalence between Assertion 
(ii), and then (i), (iii), (iv), is equivalent to the convergence of the series soon 

as forthcoming Assumption (12) is satisfied. 

Proof of Proposition 2.5. Assertion (i) implies (ii), and (ii) and (ill) are equivalent according 
to [7, Prop 7.10]. We now prove that (iv) is equivalent to (ii) and that (ii) implies (i). 

First, we check that (iv) implies that X comes down from infinity, which means that +oo is 
an entrance boundary. Then it well known that (ii) holds (see Section 8.1 in [3] or Proposition 
7.10 in [7]). Indeed, taking a = 1 in (iv), we have M := supfc>fcj E^ (exp(TfcJ) < +oo. Then, 
Markov inequality ensures that for all k > ki and t > 0, Pfc(rfc^ < t) > 1 — exp(—t)M. 
Choosing t large enough ensures that the process comes down from infinity. 

We then prove that (ii) implies (iv) by adapting the proof of [7, Prop 7.6] to the discrete 
setting. We fix a > 0 and using S < +oo, there exists ka > 1 such that 



We now define the Lyapounov function Ja as 



\ 


We notice that Ja is non-decreasing and bounded and we introduce the infinitesimal generator 
L of X, defined by 


L{f){n) = (/(n -hi) - /(n)) A„ -h (/(n - 1) - /(n)) Hn, 


for any bounded function / and any n > 1. Then, the process 
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is a martingale with respect to the natural hltration of X. Adding that LJa{m) = — 1 for any 
m > ka and that Ja{X{u)) < Ja(oo) < 1/a , we get for all k > ka and t > 0, 


Efc A TkJ)) = Efc H {aJa{X{u)) + LJa{X{u))) du) + Ja{k) 


= lEfc 




(aJa(A(n)) - 1) dn + Ja{k) 


< Ja{k). 


Adding that for any k > ka, Efc-a.s. Ja{X{t A TkJ) > Ja{ka), we get E^ 

Then (iv) follows from the monotone convergence theorem and Assumption (ii). 

It remains to show that (ii) implies (i). On the one hand, according to (2), Eoo(T/) = 
Ej_|_i(Tj) and Assumption (ii) entails that Eoo(T/) vanishes as n —>• oo as the rest of 
the hnite series S. On the other hand, under Poo, the sequence {Tn)n>o decreases to some 
random variable %oo). Then, from the monotone convergence theorem, Eoo(7ri) decreases 
to Eoo(T[o^oo)) Eoo(T[o^oo)) = 0- It ensures that T[q oo) = 0 Poo a-s. and X instantaneously 
comes down from inhnity. The proof is then complete. □ 


2.3 More tractable conditions 


Let us give some tractable conditions ensuring (1), (3), (4) or (5), which will be useful for 
examples and applications. 

Lemma 2.7. (i) Under Assumption 


k'n 

sup - < + 00 , 

n,i>l t^n+i 


limsup- < 1, 

n^oo 


( 12 ) 


Condition (3) holds if and only if 




n>l 


- < + 00 . 

hn 


(13) 


(ii) Assuming that 


hn 

sup - < + 00 , 

n,i>l Un+i 


— -^ 0, ~ ^ 

Un ^ Un 


then (1), (3), (4) and (5) are satisfied and 




kl 


n^+oo 


Un+l 


for k = 1,2, 3. 


(14) 


Criterion (13) can be seen as the discrete counterpart of the criterion in [7, p.l953] stating 
that the Feller diffusion process Z dehned by dZt = ^J^ZtdBt + Zt{r — f{Zt))dt (for a suitable 
function / and r > 0), comes down from inhnity if and only if 
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Proof. We begin with the proof of point (i). Coming back to (2), the first term of the series 
giving Kn+i{Tn) is 1/nn+i, hence E„+i(r„) > 1/nn+i- Moreover, using the second part of 
Assumption (12), there is /' < 1 such that for n large enough, \n /£ [0,^0 then for n 
large enough 

— <E„+i(ro < sup 

A^n+l l^n+j t t jJ-n+l n>l,k>0 ^J'n+k 


Then the first part of Assumption (12) allows to get (i). 


Under the assumptions (14), the properties (1), (4) and (5) are obvious, whereas (3) is a 
consequence of point (i) of the lemma. 

To get the asymptotic behavior of the moments of Tn, we use the expression of En+i(T^) 
provided in Proposition 2.2 and the fact that Xnipn goes to 0. Then, for k G {1,2,3}, by 
induction we can write 

^n+l{Tn) = k\bn,k{^ + ^n,k)i 


where hn,i = l/pn+i, bn ,2 = (E„+i(r„))2 and bn ,3 = E„+i(r„)Var„+i(r„) and An,k -t 0 as 
n —>■ 00 , which will complete the proof. Indeed, for fc = 1, we know from (2) that E„_|_i(T„) = 


1 


A^n + l 


1 + Zli 


2>n+2 


. Moreover for every I' G (0,1) and n large enough, we have 


Xnlk‘n< I' and 


An+l • • • Aj-l ^ ^/j_(n+l) 


/^n+2 ' ' ' IJ'i 


sup 


Pn 


n>l,k>0 k'n+k 


which ensures that E„_|_i(T„) ~ 1/^n+i- Combining this equivalence and the expression of 
E„_|_i(T^) provided in Proposition 2.2 yields similarly the asymptotic behavior of the second 
moment (A: = 2) and then the third moment {k = 3). □ 


Remark 2.8. Our original motivations for considering the coming down from infinity of birth 
and death processes are the regulation of large populations due to competition and the short 
time behavior of branching coalescing models (see e.g. [15] for some motivations for ancestral 
graphs). In this context, the birth rate is usually linear, whieh corresponds to independent 
reproduction events, or even zero for pure coalescing models. The death rate is often quadratic 
such as for Kingman coalescent and logistic competition, but polynomial death rate may he 
relevant, see in particular [16] for a statistical study of the death rate due to competition. 
Thus, we are interested in the particular case \n < Cn for some C > 0 and = n^log'^n 
with p > 1. In this case. Assumption (14) is obviously satisfied and then (1), (3), (4) and (5) 
hold. Proposition 2.5 ensures that the process comes down a.s. instantaneously from infinity 
and has bounded exponential moments. We refer to Section 5 for the fine description of this 
coming down from infinity. 


3 Asymptotic behavior of under Pqo 

From now on, we consider sequences {Xn)n>o and {pn)n>o satisfying the hypotheses (1) and (3). 
Thus, according to Lemma 2.1 and Proposition 2.5, Eoo is well-defined and X strongly comes 
down from infinity. Moreover T„ < -|-oo Eoo a.s. for any n > 0. In this section, we study the 
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asymptotic behavior of as n —>■ +oo under Poo- Let us recall that Eoo(7n) = '^i>n 
so that (2) yields 


^oo{Tn) = 

i>n 


1 


E 

j>i+i 


TT 


3- 




Then, S < +oo ensures that Eoo(7n) decreases to 0 as n —?• +oo. 


In the following two subsections, we compare to its mean Eqo (7n) as n —)• +oo. Two 
regimes appear depending on whether the ratio of mean times En,+i(r,i)/Eoo(7n) converges to 
a non-degenerate value or vanishes. In the first case (fast regime - Theorem 3.1), the process 
comes down very quickly from infinity, T„ is then essentially the time spent close to n and 
renormalizing T„ by its mean yields a random limit. In the second case (gradual regime - 
Theorem 3.3), can be seen as the contribution of a large number of independent random 
variables and the limit equals 1. 

In both cases, the proofs rely on the fact that Poo-a.s., where for n > 0, the 

random variable Tn is the time spent between T^+i and : 


Tn := inf{t > r„+i;X(t) = n} - Tn+i- 


By the strong Markov property, the random variables (rj)j>o are independent (under Poo) and 
Tj is distributed as Tj under P^+i. In the sequel of the section, we use for n > 0 the notation 


THn ■— E(T)2) — Efi_|_i (T)i), Tn 


_ E^^i(T'}^) 

Eoo(Tn) “ Eoo(T„) ■ 


Examples which illustrate the two regimes and the two convergences are provided in forth¬ 
coming Section 3.3, while an application to the regularly varying case is developed in Section 


5. 


3.1 The fast regime 


Theorem 3.1. We assume that (1), (3) and (4) hold and 


IEn+l(Tn) _^ 

lEoo(Tn) 


with a G (0,1]. Then, 


(d). 


lEoo(Tn) n-)-|-oo 


Z:= J]a(l-a)"Zfc, 


fc >0 


where {Zk)k is a sequence of i.i.d. random variables whose common Laplace transform G{a) := 
Eoo (exp(—a.^o)) is the unique function [0, -|-oo) —)• [0,1] that satisfies 


Va > 0, G{a) [/(l — G(a(l — a))) -|- 1 + a(l — 1{1 — a))] = 1. 


(15) 


We note that when a = 1, Z = Zq \s an exponential random variable with parameter 1. 
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Example-. If = (n!)"^ with 7 > 0, Eoo (T„) ~ ((n+1)!) Hence, lim„^+oo E(rn)/Eoo (T„) = 
1 and Theorem 3.1 (i) yields 

((n + l)!rr„ ^ E, 

n^+oo 

where E is an exponential r.v. with parameter 1. Another example is studied in forthcoming 
Section 4.2. 


Before proving Theorem 3.1, let us show the following key lemma, which focuses on the 
asymptotic behavior of the distribution of (Tn)n- 

Lemma 3.2. //lim„_,._|_oo Vn = a & (0,1], we have 

Tn (d)^ ^ 

Uln n ^+00 

where the Laplace transform off is the unique solution of (15). 

Proof. Recalling Xn/fJ-n ^ I as n —)> 00 , let us first check that 


,. Eoo(Tjj_|_i) 1 

iim ——77Z-— = lim - = 1 — a, iim n„mn-i = --77--r. 

n-s>+oo Eoo(7n) n-s>+oo mn n^+oo 1 — /(I — a) 

The first part of (16) comes from Eoo(T'„+i)/Eoo(T'„) = 1 — and 

ITT'n+l _ “^n+l Eqo (7n-|-i) 

Tn Eoo(Tn) 


(16) 


Moreover, differentiating (7) at o = 0 yields 


1 


^ _ Xn irin ^ 

h'n TUn—l hnTTln—l 


and using (4) gives the second part of (16). 

Let us prove the uniqueness of the function satisfying (15). For any bounded function g : 
[ 0 , +oo) —)■ [ 0 , 1 ], we define the function H{g) : [ 0 , + 00 ) —)• [ 0 , 1 ] as 

H{g) : a 1 —>■ [1 + a(l — Z (1 — a)) + 1(1 — g{a{l — a))]”^ • 

For two functions gi and g 2 and any a > 0, we have 


\H{gi){a) - H{g 2 ){a)\ = H{gi){a)H{g 2 ){a)l \gi{a{l - a)) - g 2 {a{l - a))| 


and using that for any a > 0, H{gi){a) < 1, 

\\H{gi)-H{g 2 )\U<l\\gi-g 2 \\oo, (17) 

which ensures the expected uniqueness since / < 1 . 
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We now prove the convergence in distribution of Tn/frin as n —)> +oo. For n > 0, let Fn : 
[0, +oo) —> [0,1] be defined as 


Fn{a) := E(exp(-ar„/mn)) = E„+i(exp(-ar„/mn)) = - , (o > 0) 


rrir 


By (7), for all o > 0 and n > 1, we have 


Gn-l 


a 


run-i 


1 + 


a 




+ -I1-G,, 


a 


rUn-l 


which we rewrite as 

Fn-l = Hn{Fn), 

where for every function / : [0, oo) —)• [0,1], n > 1 and a > 0, 


(18) 


Hnif){a) = 


1 + 


a 






m-a-l 


1 -1 


Using (16), we have for every a > 0, 


snp \Hn{f){a) - H{f){a)\ 0, 

f&Cl 

with Cl := {/ £ C^([0,oo), [0,1]) : ||/'||oo < !}• Moreover, Fn = Hn+i o ... o Hn+k{Fn+k) and 
by triangle inequality 


Fnia) - H°'‘{Fn+k){a) < F„+i(F„+i)(a) - H{Fn+i){a) 


+ 


H{Fn+i){a) - H{H^>^-\Fn+k)){a) 


Adding that for every n, Fn £ C\ and recalling (17), we get by induction over k > 0 that 


Fn{a)-H^’^{Fn+k){a) 0 

n^+oo 


for all a > 0 and A: > 0. We use again (17) to obtain that 

\\H°'^{Fn+k) - < l’^\\Fn+k - llloo < l^. 


Recalling that / < 1, we can combine the two last displays and for each e > 0, we can find k 
such that for n large enough 

|F„(a)-77°'=(l)(a)| <2e. 

Thus, {Fn{a) : n > 0) is a Cauchy sequence and Fn{a) converges to F{a) on [0,cx)). The fact 
that Fn £ cl ensures that this convergence is uniform in each compact set. Letting n —)• oo 
in (18) then yields F = H{F), which means that F satisfies (15). 

Finally, we check that F is the Laplace transform of some random variable by proving that 
F(0+) = lim(j-s.o -^(a) = 1- From (15), T(0+) is a solution of I F(0+)^ “ (1 + 0-^(0"'") + 1 = 0. 
If / = 0, this equation has the unique root 1. If ^ > 0, the two roots are 1 and 1/L But 1/1 > 1 
and obviously F{0~^) < 1, so that T(0'*“) = 1. That ends the proof of the weak convergence 
of Tn/rUn- □ 
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We can now proceed with the following proof. 


Proof of Theorem 3.1. Let Z = ~ be dehned as in the statement of the 

theorem. We need the following elementary resnlt which can be proved thanks to a simple 
indnction: for every a > 0 and for all complex nnmbers zi, Z 2 ,..., Zn,ui,.. .Un with modnlns 
less than 1 

n n 

n “ n 

i=l i=l 


<'^\zi- Ui\ . (19) 

i=l 


Then, recalling that where the r^’s are independent, 


Eoo exp -a 


Eoo(Tn 


— E (exp {—aZ)) 


E ( exp 


k>n 


—a 


Tk 


Eoo(7h) 


E f exp f —aa(l — a)^Zk 


k>0 


k>0 


E ( exp ( —a 


Tk+r 


^oo{Tn 


^ — E ^exp aa(l — a)^Zj^^ 


( 20 ) 


From Lemma 3.2, we know that in Eoo-distribntion, r^/mn converges to C- Then, thanks to 
(16) and the fact that —)• a, we have for k > 0 


'^k+n 

Eoo(Tn) 


'kkln+k TT Eoo [Tn-(-j] 
Eoo(T„+fc) Eoo[7h+i_i] 


^n+fc 



n^+oo 


a(l 


a)"C. 


The nniqneness in (15) ensnres that the variables {Zk)k are distribnted as C,. Then, with the 
last display, we get that all the terms of the snm in (20) vanish as n —>• +oo. We proceed by 
bonnded convergence. Using that 1 — exp(—x) < x for any x > 0, we get for A:, n > 0 


E ( exp ( —a 

< 

< a 


Tk+r 


Eoo(Tn) 

1 — E fexp {—a 

TYlk-\-‘ 


- E (^exp (^-aa(l - 

Eoo®)) h I' “ 


Eoo(Tn) 


+ aa(l-a)''E(Zo). 


By differentiating (15) at 0, one hnds Eoo[-^o] = 1- Moreover, 


_ Eqo (T^_|_]^) Eqo (T^_|_2) Eoo(T,i+fc) 

Eoo(r„) “ Eoo(T„) Eoo(r„+i) ■" Eoo(T„+fc_i)Eoo(r„+fc)' 


( 21 ) 


Since mfc+„/Eoo [Tfc+„] < 1 and Eoo(Tk+i)/Eoo(T„) —)-l — a<lasn—)■ + 00 , there exist 
no G N, /3 < 1 and C > 0 snch that ?nfc+„/Eoo(7n) < Cj3^ for all /c > 0, n > no. Thns, coming 
back to (21), for n > no, we have 


E ( exp ( —a 


Tk+r 


^ — E ^exp aa(l — a)^Z]^ ^ 


Eoo(Tn) 

Since the r.h.s. in the last display is snmmable, the proof is complete. 


< + aa(l-a)^. 


□ 
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3.2 The gradual regime 


We now focus on the second regime and specify the fluctuations of In this case, we will 
obtain a weak and a strong law of large numbers. 

Theorem 3.3. We assume that (1), (3) and (5) hold. 

(i) If ^ ^ 0. then 


lEoo(Tn) »^^+oo 


in Poo — probability. 


(ii) If 


and if 


Var^+i(r^) 

TL ^+oo Varoo {Tn) 


lim Varoo (T„) 

n^+oo 


y~] lEfc+idTfc — Efc+i(rfc)|^) = 0 , 

k>n 


we have 

Tn - Eqo (Tn) (d) 

Varoo (T„)V2 n^oo 

where M follows a standard normal distribution. 


( 22 ) 

(23) 


Using Cauchy-Schwarz inequality, we note that under (5), Assumption (22) implies that 
E„+i(T„)/Eoo(T„) —)■ 0, as n tends to infinity. 


Proof of Theorem 3.3. (i) - We suppose here that —)• 0. Let e > 0. Using BienaymA 

Tchebychev inequality and the independence of the random variables {Tn)n, we have 


( Tn \ Varoo(T^) ^ 

V Eoo(r„) ^ V “ e2Eoo(T„)2 e^E^iTn)^ 


(24) 


As Eoo(Tn_|_i)/m„ = l/rri — 1 —t +oo as n —)• + 00 , for all A > 0, there exists an integer no 
such that, for n > uq, Eoo(7n+i) > Anin and 

Eoo(T„)'^=(X; > 2 ^ m-fc E mi > 2A 

k>n k>n l>k k>n 


since Yli>k ~ IEoo(Tfc_|_i) > Am^. Coming back to (24), for n > no, we have 


Eoo(Tn 


- 1 


>£ < 


^ Var(rfc) 

1 k>n 


2Ae^ 


E 

k>n 


ml 


(25) 
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By assumption (5), there exists C > 0 such that Var(Tn) < Cm^. Hence, the r.h.s. of (25) 
goes to 0 as H —)> +oo and the proof of the convergence in probability is complete. 

(ii) - We follow classical ideas for the proof of central limit theorem for partial sums of inde¬ 
pendent random variables (see Theorem 27.2 in Billingsley [5]). 

We first note that Varoo (T„) > Varoo (T^+fc) for any n,k > 0 since T„ = Yli>n T the r.v. 
(ti : i > 0) are independent. Using (22), it ensures that 


Var(rfc+n) . Var(Tfc+n) 
sup —-—— < sup 


fc>o Var^o (Tji') k>o Var^o (T)i_|_fc) n—>--1-00 

This convergence being uniform with respect to /c > 0, we have 


0 . 


(26) 




1 - 


k>0 


f Var(rfc+^ 

2 (^n) J 


rsj — 


E 


Var(rfc+„) 




2 ^ Varoo (Tn 
k>o °° '' " 


Therefore, 


exp(—1^/2) = lim IT ( 1 

n^oo -1- \ 


fc >0 


Var (rfc+n) \ 

2 Varoo {Tn)) 


(27) 


'J-< _|g frp \ 

Let us now prove that converges in distribution as n —)• -|-cx) toward a standard 

normal random variable. By Levy’s theorem and (27), it suffices to prove that for any hxed t, 

f Var(Tfc+„)' 


C„ = E»(exp(V” 


HL 


fc >0 


2 Varoo {Tn) 


Varoo {Tn)^^^l 

vanishes as n —)> -|-oo. First, since the t„’s are independent, for all t G M,n > 0 

Var(rfc+ 


\Un\ = 


n TU ( ( ■,'^k+n ^{Tk+n) 


HL 


k>0 


2 Varoo {Tn) 


(28) 


According to (26), for n large enough and for any k, all the factors of the second product of 
(28) are less than 1. Hence, thanks to (19), we have the inequality 


\UA < E 

fc >0 


E fexp f 

V V Varoo (T,)^/" y, 


- 1 + 


Var(Tfc+n) 


2 Varoo {Tn) 


(29) 


According to equation (27.11) in [5, p.369], for any centered random variable ^ with a finite 
second moment, we have |E (exp(it^)) — 1-|-Var(^)t^/2| < E (min(|t^p, |t(^p)) for t > 0. 
Using this inequality with the random variables [r^+fc — E(r,p+fc)]/Varoo (T„)^'^^, we obtain 
from (29) that 

E {\Tk+n - E(rfc+n)|^) 

Varoo {Tnf^ 


\Un\<\t\^^^ 


k>0 


and using (23), Un goes to 0 as n —>■ -|-oo. This completes the proof. 


□ 


Let us now state a strong law of large numbers. 
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Theorem 3.4. We assume that (1), (3), (5) and ( 6 ) hold. Then the sequence 
converges to 1, almost surely. 


Proof of Theorem 3.4- We prove the a.s. convergence when ( 6 ) holds. According to the law 
of large numbers of Proposition 1 in [11], we just need to check that 


E 

n>0 


Var(rn) 

Eoo(7h)2 


< + 00 . 


(30) 


Using Var(Tn) < C'(E„_|_i(T„))^ thanks to (5) and Assumption ( 6 ) ensure (30) and the proof 
is complete. □ 


Let us illustrate this result with the example = 0 and = nlog'^'n (7 > 1). One can 
check that 

En+l(Tn) = 7777 77 777 

(n + 1 ) log 7 n + 1 ) 

and 

Eoo { Tji ) = ^ ^ 77 77 ~ - —7 737 . 

k>ff+l ^ ^ n^.+oo (7 _ 1 ) logT" n 

Using Lemma 2.7 (ii), we know that (1), (3), (5) hold. It’s also easy to check that ( 6 ) is true. 
Then we can apply Theorem 3.4 to get that T„/Eoo (Tn) converges a.s. to 1 as n —>■ 00 . 

Other examples will be developed in Section 5. 


3.3 Comments and examples 

In the statement of the three previous theorems appear different assumptions. Let us show 
here that the choice of these assumptions is very subtile and illustrate our results. 

1- One can exhibit a situation of gradual regime where the assumptions of Theorem 3.3 are 
satished, while (4) fails. 

We assume that for each n > 0, fin = and 

n 

Xn = Y if « e N - 4N; A„ = 2n^ if n G 4N. 


Then 


limsup An/^ri > 1 - 

n^oo 


For each n > 1, 

Y+T/2 - - 2^/2’ ^222+n/2 - - „22n/2 ’ 

SO that (1) is satished and (2) gives 


IErt+l(Tn) 
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Then (3) also holds, Ej_|_i(Tj)/E„+i(r„) is bounded for i > n, and (10) can be easily checked 
since < ^(i^^rf/ 2 - Then (5) is also fulfilled. Thus the assumptions of Theorem 3.3 are 

satisfied. 


2 - The assumptions for the weak law of large numbers in Theorem 3.3 are not sufficient to 
obtain the strong law of large numbers (Theorem 3.4). Let us consider a pure death process 
with fj,n = exp (n/log n) log n and prove that the convergence holds in probability but not 
almost surely. 

Here I = 0 and 

lE('rn) = -, Eoo (T„) = Sn+l = -. 

Moreover, as is non-decreasing. 


/ 


oo ^-x/\og{x) xoo ^-x/\og{x) g-n/log(n) 

—i-dx <Sn< —^-dx H-^- 

n log X Jn log X log n 


and 


L 


oo -x/log(x) 


-dx 

logX n->-+oo 


1 


+ 


1 


-3;/log(3;)^^ _ g-n/log(n)^ 


^logx (logx)^^ 

Combining the two last displays and recalling = E(r,i)/Eoo (Tn), we have 


Sn ~ exp(-n/logn), r„~l/logn, 0, 


so that Tn/sn+i goes to 1 in probability. 

We prove now that the almost sure convergence does not hold and proceed by contradiction. 
Thus, we assume now that Vn := T^/sn+i does converge a.s. toward 1. We have 


K+i - K = Kx+i (l- —) - —. 

\ ■Sn+l / ^n+1 

By hypothesis, the left hand side of the latter a.s. vanishes as n ^ +oo. Moreover, simple 
computations lead to Sn+i/sn —t 1 and the first term in the r.h.s. of the last display a.s. 
goes to 0 since our assumption implies that a.s. (Vn)n is bounded. Hence, putting all pieces 
together, the term Tn/sn+i has to go to 0 a.s. 

To get a contradiction thanks to Borel-Cantelli’s lemma, it suffices to prove that for e small 
enough, 

y]p(rn/s„+i > e) = oo, 

n>0 

recalling that the random variables Tn are independent. The law of Tn is exponential with 
parameter fin+i- Then, E(r„/s„+i > e) = exp(—e/r^+is^+i). Since Sn ~ logn as n —)• -|-oo, 
there exists C > 0 such that 

F{Tn/Sn+l >e)> = 4?, 

which completes the proof since X]n>o ^ (Ti-i-i/sn > e) is infinite as soon as e is small enough. 
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3 - In Theorems 3.1 and 3.3, we did not consider the case where = E(r„)/Eoo (Tn) does 
not converge. In such a case, one can only state analogous results along the convergent 
subsequences. For instance, if fj, 2 n = fJ- 2 n+i = 3^", we have 


rori —- and 

n^+oo 9 


r2n+l 


n^+oo 


4 

5' 


Theorem 3.1 then still holds but the subsequences (T2n/Eoo {T 2 n))n and (T2„+i/Eoo {T 2 n+i))n 
converge in distribution to different limits. 

One can also find examples where 0 = liminf^r^ < limsup„r„. Then, (T„)„ has two subse¬ 
quences satisfying the two regimes of Theorems 3.1 and 3.3. 


4 Speed of coming down from infinity 

In this section, we use the asymptotic behavior of T„/Eoo(T'„) obtained in the previous Section 
to derive the short time behavior of X[t). We prove that X behaves as the following non¬ 
increasing function that tends to infinity as t —>■ 0 

v{t) := inf{n > 0; Eoo(T„) < t}- 

The function u is a cad-lag step function defined on M_|_, decreasing from -|-oo to 0, it equals 
n between Eoo(T„) and Eoo(T„_i), and 0 after Eoo(To). 

The short time behavior of X relies on the inversion of the asymptotic behavior of T„ (see 
forthcoming Lemma 4.2) and the control of the excursion of X between two successive stopping 
times [r„+i,T„] (see forthcoming Lemma 4.1). The latter is true under the assumption 

limsup — < 1. (31) 

n^oo 

This assumption is already necessary for Theorem 3.1 but not for Theorem 3.3, as developed 
in Subsection 3.3 Example 1. 

The proof is organized as follows. We introduce the a.s. non-increasing process Y defined by 

Y{t) = n if te[Tn,Tn-i) 

In the next Section, we prove that this (more regular) process comes down from infinity at 
speed v(t) and we compare the processes X{t) and Y{t) as t —)• 0 by the study of the height 
of the excursions of the process X. 

4.1 Height of the excursions and non-increasing process. 

We first compare the processes X (t) and Y (t) by estimating the number of birth events 
between the times Tn and r„_i: 

Hn = #{s G [T„,T„_i) : X{s) -X(s-) >0}, n > 1. 
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Lemma 4.1. (i) We have 


X{t) Y{t) Hyjt) Y{t) 

~ v{t) v{t) ~ Y{t) v{t) 


(32) 


(a) Under Assumption (31), —- —)■ 0 Poo O'-S- 

n 

Proof, (i) For any t G [T„,T„_i), Y{t) = n and 0 < X{t) — Y{t) < Hn, so the first part is 
obvious. 

(ii) Let us first notice that equals the number of positive jumps between time and 
r„_i of a random walk whose transition probabilities are given by Pi^i+i = + tii), 

Pi^i-i = /ii/(Ai + ^i) for i > 1. Using (31), we can choose no large enough so that p = 
sup„>^Q Xn/{^n + hn) < 1/2- Then, for n > no, Hn is stochastically dominated by T, the 
hitting time of n — 1 by a simple random walk starting at n, with probability transitions 
(1 —p,p). Since p < 1/2, E(T^) < +oo. Hence sup„>„g Eoo (^n) < +oo and the sequences 
(Eoo iHn))n and (Eoo {Hn))n are bounded. 

Let us now consider the Laplace transform of Hn given by Gnia) = Eqo (exp(—aFf^))- la the 
same vein as we have obtained (7) and by applying the strong Markov property at the first 
time when X jumps after T„, we get the recursion formula 

Gn{n) = -—— -h -— -e °'Gn{ci)Gn+i{ci), a P 0, n > 1. (33) 

An Y hn An + hn 

Differentiating (33) twice at o = 0, the second moment of Hn satisfies the following recursion 
formula 


^Eoo {Hi) = Eoo (FF^+i) + 1 + 2 (Eoo {Hn) + Eoo (FFn+l) + Eoo {Hn) Eoo (FFn+l)) • 

We have seen that the right hand side of the latter is uniformly bounded in n > 0. It entails 
that there is C > 0 such that 

Eoo(FF 2) <c^, n>l. (34) 

f^n 

Finally, Eoo (j2n>i ^ ^ Y.n>i (31)- la particular, it turns 

out that the sequence (^)n almost surely goes to 0 as n —?• +00. □ 


Let us now introduce the quantity 

R{x,y) := Eoo(r[,j,])/Eoo(r[y]) 
and study the behavior of Y{t)/v{t) as t tends to 0. 

Proposition 4.2. (i) //Tn/Eoo(Tn) is tight on (0, 00) and for every x > 1, limn->.oo R{nx, n) = 
0 then Y{t)/v{t) —>• 1 m Poo — probability. 

(ii) //rn/Eoo(Tn) 1 in Poo — probability and if for every x > 1, limsup^_,.oo F?(nx, n) < 1, 
then Y{t)/v{t) —?• 1 m Poo — probability. 
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(iii) //lim„_^oo R{n + l,n) = 1 and limsup„_^oo R{nx,n) < 1 for every x > 1, then 

lim limsup I _ l| = Q. (35) 

a^l 4^0 v{t) 

If additionally r„/Eoo(7n) —>■ 1 a.s., then Y{t)/v{t) —)• 1 Poo cl-S- 

Proof, (i) Under the tightness assnmption, for any e > 0, there exist 0 < ^4 < i? snch that for 
every n > 0, 

Poo(^<r„/Eoo(r„) <S) > 1-e. (36) 

Moreover, for every x > 1 and n large enongh, 

R{nx,n) <mm{l/{2B), A/2). (37) 

By the dehnition of the fnnction v, we have 


lEoo {Tvit)) <t <Eoo (T^(i)_i) 


It implies that for any t > 0, 


'oo ( T[^v(t)] - 2 ) - 1 - 


Eo 




T[xv{t)] ^ R{v{t),xv{t)) 


j ^Eoo {T[^vit)]) 


Hence, nsing (36) and (37) and for t small enongh. 


We similarly get that for t small enongh 


< B\ > 1 - e. 


Poo(r[,p)/,,] > 2t) > Poo I ^ s > AI > 1 - e. 

Then, we have for t small enongh 

IF’oo(T'[j,„(i)] < t/2,r[^p)/3,] > 2t) > 1 - 2e. 


Since Y is non-increasing, that implies Poo(U(t) G [v{t)/x,v{t)x]) > 1 — 2e and ensnres that 
Y{t)/v{t) tends to 1 in probability as t —)• 0. 

The proof of (ii) follows the same steps as the one of (i). Since rn/Eoo(7n) —)• 1 in Poo — 
probability as n tends to inhnity, we can choose for any x > 1, A and B close enongh to 1 
and a < 1 snch that 


Poo ( ^ < - Kb/^] —^ H I > 1 — e, R{nx, n) < a min(l/H, A) 

\ IEcxd [T[v(t)/x]) J 

for t small enongh and n large enongh. We conclnde as previonsly nsing now Poo(T'[3,^(^)] < 
^ij'R[v(t)/x] — — 1 2e. 
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Let us now prove (iii). We first note that v is non-increasing. Let us first prove that 


limliminf = i 

a 4 ,l t^o v(t) 


Then for all t > 0 and a > 1, 


yjat) 

vit) - ■ 


( 38 ) 


Moreover, 


IEoo(r^(ai)) ^ IEoo(T^(i)) R{v{t) - 1, v{t)) = Eoo(T„(i)-i) > t. 

For all r] > 1, the hrst assumption of (iii) yields R{v{t) — l,v{t)) < rj for t small enough and 
we get 

IEoo(T;;(aL) ^ «hIEoo(r^(t)), 

which implies that 

limlimsupii(n(at),n(t)) < 1. 

Moreover the second assumption of (iii) ensures that liminf^^oo "ra) > 1 for every 

X < 1. We add that for each t such that v{at) < xv{t), we have 

R{v{at),v{t)) > R{xv{t),v{t)) 

Combining the three last displays ensures that for any x < 1 ensures that 

limliminf > x. 

ail t-s-0 v{t) 


Letting x —)• 1 and recalling (38) yields 


,v(at) I 

hm hm sup —-1 = 0. 

a;i ' v{t) 


To conclude to the hrst part of (iii), it remains to consider a < 1, which is simply derived 
from the previous limit by changing t into t/a. 

We assume now that r„/Eoo(Tn) —>• 1 a.s. and (35) enables us to compose the equivalence by 
v: 


since u(Eoo(T„)) 
to 0 a.s. and 


vjTn) 

u(Eoo(r„)) 


= 1 = lim 
n^oc 


v{Tn 


a.s. 


n 


n by dehnition of v and n —)> Eoo(2~n) is decreasing. Noting that Tn goes 


n y(t) n 

v{Tn) ~ v{t) ~ v{Tn-i) 

a.s. on the event t G [r„,T„_i) ends up the proof of (iii). □ 


21 











4.2 Fast coming down from infinity 

We state the convergence in probability inherited from Theorem 3.1 (a > 0). 
Theorem 4.3. We assume that (1), (3) and (4) hold and 


En+i{Tn)/Eoo{Tn) aG(0,l]. 


Then 

X(t) 

lim — = 1 in Poo — probability, 
t^o v{t) ^ " 

Proof. Under the assnmptions of Theorem 3.1, the seqnence r„/Eoo(7n) converges in law to 
a random variable whose law is snpported on (0, oo), since G{a) —)• 0 as a —?• oo. Then this 
seqnence is tight on (0, oo). Let ns fix x > 1 and show that 


lim R{nx, n) = lim 

n^oo n—¥oo 


IEoo(T[„a;]) 

Eoo(Tn) 


= 0 . 


Indeed, as seen in (16), lim„^_|_oo 


= 1 — 0. Then for n large enongh. 


lim = lin, 'tt'Eoo(^ 

n^oo Eoo(T„) n^oo ii Eoo(Tj) 


lim (1 -a) [”^1-^ = 0. 

n^oo 


Then by Proposition 4.2-(i), Y{t)/v{t) —?■ 1 in Poo — probability. Lemma 4.1 will then allow 
to conclnde. □ 


Example. Let ns consider a pnre death process with p.n = with ,0 > 0 and thns Eoo (7n) ~ 
!{\ — e“^). Hence, the conditions of Theorem 3.1 hold trne with a = 1 — and 


n^+oo 


Y^e-^'^Ek 


where the E}fs are independant exponential random variables with parameter 1. In that case, 
we can explicitly determine the speed v of Theorem 4.3 and we get X{t) -{\ogt)/j3 as t 0 
in probability. 


4.3 Gradual coming down from infinity 

We give now the speed of convergence and describe the flnctnations of X in the case a = 0. 
Theorem 4.4. ITe assume that (1), (3), (5) hold and E„_|_i(r„)/Eoo(Tn) —> 0- 

n^+oo 

ITe assume also that (31) holds and that for every x > 1, 

lim sup < 1. (39) 

n^oo^ Eoo(r„) ^ ’ 

Then, 

X(t) 

lim —^ = 1 in Poo — probability. 
t-s-o v{t) 

Assuming further that 'Yfn (En+i{Tn)/Eoo{Tn)^ < oo, this convergence holds a.s. 
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The additional assumption (39) is required for the inversion (Proposition 4.2 (ii)). Indeed the 
quantity Eoo(Tn) should not tend too slowly to 0. The example Hn = nlog'^n (7 > 1) and 
An = 0 shows that (39) may fail while the other assumptions hold (see Section 3.2 for details). 


Proof. Let us deal with the convergence in probability and work under Pqo- The first four 
assumptions allow us to apply Theorem 3.3, so that Tn/Eoo(Tn) —>• 1 in probability as n —)• 
00. Using (39), we can apply Proposition 4.2 (ii) to get that Y{t)/v{t) —?■ 1 in probability. 
Moreover Assumption (31) enables us to use Lemma 4.1. It ensures that HY(t)/Y{t) —)• 0 a.s. 
since Y(t) —>■ +00 a.s as t —)• 0 and X{t)/v{t) — Y(t)/v(t) —)• 0 in probability and then the 
convergence in probability of X{t)/v{t) to 1. 


We note now that 


R{n + 1, n) = 1 — 


IEn+l(Tn) 

lEoo(Tn) ■ 


Thus, E„_|_i(T„)/Eoo(Tn) —?• 0 yields R{n + 1, n) —?• 1 as n —?• 00. Then the a.s. convergence is 
obtained similarly combining Theorem 3.4, Lemma 4.1 and Proposition 4.2 (hi). □ 


Remark 4.5. We remark that if Xn = 0 and fj,n = n{n — l)/2, X{t) is the number of bloeks 
of the Kingman eoalescent at time t. In this ease, 


Eoo(Tn) = 

i>n 


2 

iii + 1) 


2 

n 


Then, v{t) = j and we reeover from Theorem 4-4 ^^6 speed of coming down from infinity for 

this process, obtained by Aldous in [ll-paragraph 4-2.: tX{t) —2 a.s.. 

t—>-0 


We refer to the next section for more general examples, where we also provide the fluctuations 
of X under Pqo for t close to 0 using the following result. 


Proposition 4.6. We assume that (1), (3) and (5) hold and that (22) and (23) hold. 
We also assume that n ^ ^ every x G M, 


t lEoo {Ts{x,t)) _^ 

^Varoo (r,(,,i)) 


where s{x,t) = [v{t) + x^/v{t)]. Then, 


VW) 




t-s-O 


AA, 


where J\f follows a standard normal distribution. 


(40) 


(41) 


Proof of Proposition 4-6. The proof follows the same steps than the previous theorem. We use 
the C.L.T theorem for Tn to firstly establish a central limit theorem for the a.s. non-increasing 
process Y. 
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As V is non-increasing, we can follow the proof of the central limit theorem for renewal 
processes (as suggested by Aldous for Kingman’s coalescent, cf. [1]). More precisely, for any 
t > 0 ,x G M, we use Poo(K(t) > s(x,t)) = Poo(T's(x,t) > to get 


Poof 


m 

v{t) 


- 1 > X = 


■‘s{x,t) 


> 


t Poo (Ts{x,t)) 

^Varoo {Ts(^x,t)) j 


where 

~ _ r„ - Eoo {Tn) 

^Tl — /- 

Y^Varoo (Tn) 

All the assumptions of Theorem 3.3 (ii) are met, so Zn converges weakly to a standard normal 
variable. Using (40), we obtain the C.L.T. (41) for Y. 

We end the proof by deducing the C.L.T for X thanks to the decomposition 







X{t) - Y{t) 

VW) 


(42) 


From (32), we almost surely have 


VW) ~ \/W) 

Using (34), there exists C such that 



(43) 


Since this series converges by hypothesis, Hn/y/n a.s. goes to 0 as n —)> -|-oo. Recalling that (5) 
and (22) ensure that E„+i(r„)/Eoo(T„) —>■ 0, all the assumptions of Theorem 4.4 are fulhlled, 
so Y{t) r\j X{t) ~ r;(t) as t —)• 0 in probability. Then the right hand side of (43) vanishes 
as t — 7 - 0 in probability and (42) allows us to derive (41) from the C.L.T for Y established 
above. □ 


5 Application for regularly varying death rates 

As mentioned before (see in particular the end of Section 2.3), the following class of birth 
and death processes is particulary relevant for population dynamics and population genetics 
models: 


An < Cn, Hn = log"^ n, where C > 0, p > 1 ,7 G M. 


(44) 


This is a particular case of a main class that we can attain with our results. In what follows, 
we will suppose that the birth rate is sub-linear as assumed in (44). The main assumption is 
that the death rate varies regularly. Our previous theorems apply in this general context. 


24 




















Recall that a sequence of real non-zero numbers {un)n>o varies regularly with index /? 7^ 0 if 
for all a > 0, 


lim 

n—>-+oo Un 




A function 5 : [0, -|-cx)) 


(0, -|-cx)) varies regularly at 0 with index p 7^ 0 if for all a > 0, 
g{ax) 


lim 

x^o g[x) 


= aP. 


Theorem 5.1. Suppose that lim„^+oo = 0 and that {iJ.n)n varies regularly with index 

p > 1. Then, 

lini - —Tn = lim = 1 Poo - as- (45) 

n^oo n t-S’O v{t) 

where v is regularly varying at 0 with index 1/(1 — /o). Further, 

_^ M, 

i/n V " /r„+i(p - 1)/ n^+00 

Assuming further that X]n>i ^ +00, we also get that 

converges in law, as t tends to 0, to a standard normal distribution. 

We recover the central limit theorem for the Kingman coalescent. We also mention that [14] 
provides Gaussian limits for more general A coalescent processes whose “Kingman part” is non 
trivial. 

Corollary 5.2. Assume (44). Then, 

lim Tn {p — l)nP~^ log"^ n = 1 Poo ~ “-S- 

n^oo 

and 

_i1 Af 

y/n V ” [p — l)nP ^log'^n/ n^+ca 

where M follows a standard normal distribution. 


Proof of Theorem 5.1. Let us hrst remark that since {pn)n varies regularly with index p > 1, 
then (l//in)n (resp. {1/pf^)n) varies regularly with index —p < —1 (resp. —2p < —2). Under 
the assumptions of Theorem 5.1, we note that conditions (14) are satished. Indeed, Theorem 
1.5.3 in [6] shows that the sequence (pn) is equivalent to a non-decreasing sequence and Lemma 
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A.l applies directly to Lemma 2.7 (ii) can be applied, so that (1), (3), (4) and (5) 

are satisfied and 


^n+iirH 


r\j 

n^+oo 


kl 

h ■ 

<+i 


for fc = 1, 2, 3. 


(46) 


Then by forthcoming Lemma A.l, 


Eoo(Tn) ~ 

n^+oo 


n 

k'n+liP ~ 1) 


Thns for x > 1, 

n^oo ^oo\-^n) 

Since 1 — p < 0, then < 1. Moreover, (En+i(T,i)/Eoo(T„))^ converges and the 

assnmptions for Theorem 3.4 and Theorem 4.4 are satished, implying (45). 

To prove the C.L.T. for (T^), we nse again (46) for k = 1,2 to get Var(rn) r\j 1 /pI+i as 
n —?• + 00 , which implies that (Var(rn))n varies regularly with index —2p. Then Varoo(7n) = 
X)^ri+i and forthcoming Lemma A.l ensures that Varoo (Tk) varies regularly with 

index 1 — 2p and 


Varoo (Tn) ~ 


n 

~ ^)Pn+l 


(47) 


Therefore we have . ~ 2 p_L entails (22). 

Var«)(T„) ^ ^ ' 

Moreover, by the triangle inequality and the binomial theorem, we have 

IEfe+i(|7fc — Efc+i(Tfc)|^) < Efc+i(r|) + 3Efc+i(rfc)Efc+i(r|) + 4Efc+i(rfc)^. 


Thanks to (46), all the terms of the r.h.s. are of order of magnitude 1/p^j^i as n ^ +oo. 
Thus, using again Lemma A.l and (47), there is a positive constant C such that 

^^k>n lEfc+idTfc - Efc+i(rfc)|^) ^ C' 

Varoo (T„)"/' “ 

This latter vanishes as n —)• +oo and (23) is satished. Hence we apply Theorem 3.3 and 


_ ^aa, 

y/n V ” p„+i(p - 1)/ n^+oo 
where Af follows a standard normal distribution. 

Let us now prove the last assertion of Theorem 5.1. To apply Proposition 4.6, we need to 
prove that 

^ (48) 

y'Varoo (r,(,,t)) 

and we hrst consider the case x > 0. We note that 


lEoo {'A'v{t)) lEoo {Ts(x,t)) ^ t Eoo {Ts(x,t)) ^ {Tv(t) — l) 1^00 (2^s(a;,i)) ) 


(49) 
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and we handle the two sides similarly. For the left hand side, we have 


[n+xy^ —1 _ 

lEoo (In) — lEoo f^[n+Xv^) ~ ^k+l{Tk) — , (50) 

\ / n—>+oo 

k=n 


using forthcoming Lemma A.3 with Un = [n + x^/n]. Moreover, applying forthcoming Lemma 
A.2 with f{y) = y, g{y) = [y + u^] and h{n) = Vaioo (F„) and using (47), we get 


Vaioo 


- [n+xy^] 


~ Vaioo {Tn) 

n^+oo 


n 

rsj -^ 

n^+oo (2p- 


Combining this equivalence with (50) and (49) yields (48) for x > 0, while the case x < 0 can 
be handled similarly. It ends up the proof. □ 


Proof of Corollary 5.2. We easily remark that 


IEn+l(7n) 


1 

rsj - 

n^+cxD nP log"^ n 


Eoo {Tn 


r\j 

n^+oo 


E 

k>n+l 


kP log"’' k n->+oo (p — l)n^-i log"^ 


n 


and ^ 

“ ^ n^+oo ((2p — 1)77,2/’“’^ log^’’'n) 

The assumptions of Theorem 5.1 can then be easily checked to get the result. □ 


A Appendix : regularly varying functions 

The proofs of the previous section rely on the following technical results on regularly varying 
functions. 


Lemma A.l. Let g be a function that varies regularly at +oo with index p' < —1. Then the 
series converges and R{n) = Ylk>n varies regularly with index p' + ^ and 


'^9{k) 

k>n 


r\j 

n^+oo 


ngjn) 

p' + T 


Proof. First, since p' < —1, X]fe>o5(^) are both convergent. Moreover, 

thanks to [6, Thm 1.5.3], a regularly varying function with negative index is equivalent to 
a non-increasing function. Then, without loss of generality, one can suppose that g is non¬ 
increasing. Then, if In ■= g{x)dx, a classical comparison between series and integrals 
entails that 1 < ^ < 1 -|- Using that g varies regularly and according to [6, Thm 1.5.11], 

. nqin) , , ^ 

lim —-— = -{p + 1). (51) 

n^+oo In 

Hence, In ~ Rn ss n ^ -|-oo. We also see from (51) that I varies regularly at -|-oo with index 
p' -|- 1. Since I and R are equivalent, R also varies regularly with the same index. □ 
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Lemma A.2. Let xq G [0,+oo] and let f and g be two positive functions such that 


fix) 


>*0 


L G {0, +oo}, 


fix) , 

g{x) x->-xo 


If h varies regularly at L, then 

Hfjx)) ^ ^ 

h{g{x)) x^xo 

Moreover, if fix) = /(x,t) = gix){l + te{x)) with lima;_>.a;(j e{x) = 0, the previous convergence 
holds uniformly in t in any compact subset o/M. 


Proof. We only prove the case L = 0 and fix e > 0. Thanks to Theorem 1.5.1 p.22 in 
[6], the convergence given above in the definition of regnlarly varying fnnction can be taken 
nniform with respect to a in some compact set. Then, there exist r/' > 0 snch that for every 
a £ [I — r],l + g] and y G (0, g'), 


l-e< 


hjay) 

Hy) 


< 1 + e. 


Fnrthermore, for x close enongh to xq, we have g{x) < rf and (1 — ?]) < fix)jg{x) < (1 + g) 
so that 



< e, 


which ends np the first part of the proof. The second part follows in the same way since 
1 + te{x) goes to 1 nniformly in t in any compact set. □ 


Lemma A.3. Let imn)n be a regularly varying sequence and iun)n a sequence of integers such 
that Un ^ oo and Un/n —>-0 as n —>■ oo. Then 


n+Un—l 

ruk 

k=n 


~ UnTUn. 
n^+co 


Proof. We write 


ri+Un-l 


ruk 


- Unmr 


k=n 


n+Un-1 




k=n 


Ulk 


rUr, 


- 1 


< UnrUn sup 
te[o,i] 




nir. 


Using the second part of Lemma A.2 with the regularly varying sequence imn)n, 
vanishes, which ends up the proof. 


- 1 


the last term 
□ 


Acknowledgement. This work was partially funded by the Chaire Modelisation Mathema- 
tique et Biodiversite VEOLIA-Ecole Polytechnique-MNHN-F.X., by the labex LMH through 
the grant no ANR-11-LABX-0056-LMH in the "Programme des Investissements d’Avenir", by 
the professorial chair Jean Marjoulet and by the project MANEGE ‘Modeles Aleatoires en 
Ecologie, Genetique et Evolution’ ANR-09-BLAN-0215. 


28 















References 

[1] D. J. Aldous. Deterministic and stochastic models for coalescence (aggregation and co¬ 
agulation): a review of the mean-field theory for probabilists. Bernoulli, 5(l):3-48, 1999. 

[2] L. J. S. Allen. An introduction to stochastic processes with applications to biology. CRC 
Press, Boca Raton, FL, second edition, 2011. 

[3] W. J. Anderson. Continuous-time Markov chains. Springer Series in Statistics: Proba¬ 
bility and its Applications. Springer-Verlag, New York, 1991. An applications-oriented 
approach. 

[4] J. Berestycki, N. Berestycki, and V. Limic. The Lambda-coalescent speed of coming down 
from inhnity. Ann. Probab., 38(l):207-233, 2010. 

[5] P. Billingsley. Probability and measure. Wiley Series in Probability and Mathematical 
Statistics: Probability and Mathematical Statistics. John Wiley & Sons Inc., New York, 
second edition, 1986. 

[6] N. H. Bingham, C. M. Goldie, and J. L. Teugels. Regular variation, volume 27 of Ency¬ 
clopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 
1989. 

[7] P. Cattiaux, P. Collet, A. Lambert, S. Martinez, S. Meleard, and J. San Martin. Quasi¬ 
stationary distributions and diffusion models in population dynamics. Ann. Probab., 
37(5):1926~1969, 2009. 

[8] P. Donnelly. Weak convergence to a Markov chain with an entrance boundary: ancestral 
processes in population genetics. Ann. Probab., 19(3):1102-1117, 1991. 

[9] S. Karlin and J. L. McGregor. The differential equations of birth-and-death processes, 
and the Stieltjes moment problem. Trans. Amer. Math. Soc., 85:489-546, 1957. 

[10] S. Karlin and H. M. Taylor. A first course in stochastic processes. Academic Press [A 
subsidiary of Harcourt Brace Jovanovich, Publishers], New York-London, second edition, 
1975. 

[11] O. I. Klesov. The rate of convergence of series of random variables. Ukrain. Mat. Zh., 
35(3):309-314, 1983. 

[12] M. Kot. Elements of mathematical ecology. Cambridge University Press, Cambridge, 

2001. 

[13] A. Lambert. The branching process with logistic growth. Ann. Appl. Probab., 15(2):1506- 
1535, 2005. 

[14] V. Limic, A. Talarczyk (2014). Diffusion limits for mixed with Kingman coalescents at 
small times. Avialable via http://arxiv.org/abs/1409.6200. 

[15] P. F. Slade, J. Wakeley. The Structured Ancestral Selection Graph and the Many-Demes 
Limit. Genetics. 2005 169(2): 1117-1131. 


29 



[16] R. M. Sibly, D. Barker, M. C. Denham, J. Hone, and M. Pagel. On the regulation of 
populations of mammals, birds, fish, and insects. Science, 309(5734):607-610, 2005. 

[17] H. M. Taylor and S. Karlin. An introduction to stochastic modeling. Academic Press, 
Inc., San Diego, CA, third edition, 1998. 

[18] J. van den Broek and H. Heesterbeek. Nonhomogeneous birth and death models for 
epidemic outbreak data. Bio statistics, 8(2):453-467, 2007. 

[19] E. A. van Doom. Quasi-stationary distributions and convergence to quasi-stationarity of 
birth-death processes. Adv. in Appl. Probab., 23(4):683-'700, 1991. 


30 



